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The kinetics of relaxation in quantum systems induced by anomalously slowly fluctuating noise is 
studied in detail. In this study two processes are considered, as examples: (1) relaxation in two-level 
system (TLS) caused by external noise with slowly decaying correlation function P(t) ~ (wt)~ a , 
where < a < 1, and (2) anomalous-diffusion controlled radical pair (RP) recombination in which 
relaxation results from the reaction with slowly fluctuating reaction rate whose fluctuations are 
governed by subdiffusive relative motion of radicals in a potential well. Analysis of these two 
processes is made within continuous time random walk approach (CTRWA). Rigorous CTRWA- 
treatment of the processes under study results in the non-Markovian stochastic Liouville equation 
(SLE) for the density matrix of the systems. This SLE predicts important specific features of 
relaxation kinetics of quantum systems in the presence of the above mentioned anomalous noise. In 
TLS, for example, relaxation of both phase and population turns out to be anomalously slow and 
strongly non-exponential. Moreover, for a < 1 in the limit of characteristic fluctuation rates w much 
larger than the frequency of quantum transitions ui a (w/uj s S> 1) the relaxation kinetics in TLS is 
independent of w. Strong changes of the relaxation kinetics is found with the change of to and a. 
In RP recombination the anomalous fluctuations of reaction rate (due to anomalous diffusion) show 
themselves in non-analytical dependence of the reaction yield Y r on reactivity and parameters of 
the RP spin Hamiltonian. In particular, the spectrum shape of reaction yield detected magnetic 
resonance, i.e. the dependence of Y r on the frequency ui of resonance microwave field is found to be 
strongly non-Lorenzian with the width to be determined by the strength wi of this field. 

PACS numbers: PACS numbers: 05.40.Fb, 02.50.-r, 76.20.+q 



I. INTRODUCTION 

The noise induced relaxation in quantum systems is 
very important process, which is investigated in vari- 
ous brunches of physics and chemistry: magnetic res- 
onance 0, quantum optics, nonlinear spectroscopy 0, 
solid state physics 0)3, etc. 

Some of these processes are analyzed assuming fast 
decay of correlation functions of this noise r c and con- 
sidering the short correlation time limit (SCTL). This 
approach is well known in the general relaxation theory 
[5| . In the absence of memory the relaxation is described 
by very popular Bloch-type equations. Nevertheless, in 
a large number of other processes the memory effects 
come to play. Important steps beyond the conventional 
short correlation time approximation can be made within 
the approach treating the relaxation kinetics with the 
stochastic Liouville equation (SLE) [|J. This approach, 
however, also has strong limitations: it assumes the noise 
fluctuations to be Markovian. The alternative approach, 
which can be applied for the analysis memory effects is 
based on the Zwanzig projection operator formalism 
Unfortunately, in reality the projection formalism allows 
for tractable analysis of the memory effects only in the 
lowest orders in the fluctuating interaction V inducing 
relaxation @. 

At present especially significant interest to the correct 
interpretation of long memory (non-Markovian) effects 
has been excited by recent works concerning processes 
governed by noises whose correlation functions P(t) de- 



cay anomalously slowly: P(t) ~ t~ a with a < 1. A 
large number of different phenomena, in which these pro- 
cesses play important role, are considered in the review 
article Some of them are discussed in relation to 
spectroscopic studies of quantum dots |nj. Sim- 
ilar problems are analyzed in the theory of stochastic 
resonances 0]. Many works study the manifestation 
of non-Markovian long memory effects in dielectric re- 
laxation and linear dielectric response |13|| (and refer- 
ences therein). All anomalous relaxation processes men- 
tioned above cannot be properly described by any meth- 
ods based on the short correlation time approximation. 
The conventional SLE approach is not appropriate for 
description of these processes either. 

The efficient method of analyzing the memory effects 
(including the case of anomalously long memory) has 
recently been developed in ref. 0|. It treats the 
noise within the continuous time random walk approach 

(ctrwa) m mill, describing the statistical proper- 
ties of the noise in terms of the probability distribution 
functions (PDFs) of renewals of the interaction fluctua- 
tions associated with the noise. The developed CTRWA- 
based method is based on the Markovian representation 
of the any CTRW-processes by the Markovian ones pro- 
posed in ref. This representation allows for obtain- 
ing the formal matrix expression for the evolution oper- 
ator of the system under study, which can be considered 
as the non-Markovian generalization of the conventional 
SLE 0] and is called hereafter the non-Markovian SLE. 

The non-Markovian SLE appeared to be very fruit- 
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ful for description of some classical processes assisted by 
stochastic anomalous spatial migration |18|. The mem- 
ory effects are known to manifest themselves in this 
processes extremely strongly leading to stro ngly non- 
exponential anomalous relaxation kinetics |l9l l20l |21| . 

In this work the proposed non-Markovian SLE is ap- 
plied to the detailed analysis of specific features of relax- 
ation in quantum systems induced by anomalous noise, 
i.e. the noise with anomalously long-tailed correlations 
in time, whose effect on quantum systems can hardly 
be correctly treated by expansion in powers of the noise 
amplitude. This type of the SLE enables one to de- 
scribe relaxation kinetics without expansion in the fluc- 
tuating interaction. In some physically reasonable ap- 
proaches, for example in the sudden relaxation model 
(SRM), it allows for describing the phase and population 
relaxation kinetics in analytical form. In the work two 
types of processes are considered: relaxation in two-level 
systems (TLS) and anomalous-diffusion assisted radical 
pair (RP) recombination in a potential well, as an ex- 
ample of processes in multilevel systems. Relaxation in 
these systems is shown to be strongly non-exponential. 
In addition, in RP recombination anomalous properties 
of relaxation (caused by anomalous relative diffusion ) re - 
sult in some peculiarities of magnetic field effects [22J: 
non-analytical dependence of magnetic field affected re- 
combination yield (MARY) on the parameters of the spin 
Hamiltonian, strongly non-Lorenzian shape of lines of re- 
action yield detected magnetic resonance (RYDMR) [2i| , 
etc. 



II. GENERAL FORMULATION 

We consider noise induced relaxation in the quantum 
system whose evolution is governed by the hamiltonian 

H(t)=H s + V{t), (2.1) 

where H s is the term independent of time and V(t) is the 
fluctuating interaction, which models the noise. The evo- 
lution is described by the density matrix p(t) satisfying 
the Liouville equation (H = 1) 

p = -iH(t)p, with Hp = [H, p] = [Hp - pH}. (2.2) 

^(^-fluctuations are assumed to be symmetric ((V) = 
0) and result from stochastic jumps between the states 
\x v ) in the (discrete or continuum) space {x u } = {x} 
with different V — V v and H = H v (i.e. different V = 
V v = [K, ■ • ■] and # = #,): 

V = V \x„)Vu{x v \ and7i; = y \x u )H v {x„\. (2.3) 

Hereafter we will apply the bra-ket notation: 

\k), \kk') = \k)(k'\, and \x) (2.4) 

for the eigenstates of H (in the original space) and H (in 
the Liouville space), and for states in {xj-space, respec- 
tively. 



It is worth noting that within the proposed semiclas- 
sical approach one can also describe some fluctuating ir- 
reversible processes in the system modelled by the addi- 
tional non-Hermitian relaxation term —iK{t) in the Li- 
oville operator H. In this model the effect of fluctuating 
irreversible process results in modification of the Liou- 
ville equation 1)2.2(1 : 

p = —iLp, in which L = TL — ilC. (2-5) 

In the majority of processes the time evolution of ob- 
servables under study is determined by the evolution op- 
erator TZ(t) in the Liouville space averaged over V (t)- and 
-ftf(£)-fluctuations: 

pit) = ti(t) Pi , where K(t) = (re-^o dri ^ . (2.6) 

The operator lZ(t) can equivalently be represented in 
terms of the conditional evolution operator Q(x, x'\t) av- 
eraged over the initial distribution Pi{x): 

n(t) = (g)=y] $(x,Xi\t)Pi{xi). (2.7) 

* — 4 X,Xi 

For steady state ^(t)-fluctuations the averaging should 
be made over the equilibrium distribution P e (x), i.e. 
Pi(x)=P e (x). 

It is worth noting, however, that evaluation of some 
observables requires analysis of the conditional operator 
Q(x,Xi\t) itself rather than the averaged one it(t) (see 
Sec. IVC). 

Thus, in the proposed semiclassical approximation for 
V(t)- and i^(i)-fluctuations the relaxation kinetics is de- 
termined by the operator Q(x,Xi\t). Its evaluation is, in 
general, a complicated problem. In what follows we will 
discuss important approaches in which this problem can 
be significantly simplified by reducing it to solving the 
(differential or integral) SLE for G(x,Xi\t) [6j. 

In our further study we will often use the Laplace 
transformation of functions under study in time t con- 
ventionally denoted as 

Z(e) = / dtZ{t)e- et (2.8) 
Jo 

for any function Z(t). 



III. STOCHASTIC LIOUVILLE EQUATIONS 

The goal of our work is to analyze the strong mem- 
ory effects within the CTRWA. This approach is known 
to result in the complicated non-Markovian SLE ^4| f° r 
Q(x,Xi\t). However, it is instructive to start the discus- 
sion with the simpler Markovian models reducing to the 
conventional semiclassical SLE. 
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A. Markovian fluctuations 

In the Markovian approach the kinetics of jumps in 
{x}-space leading to ^(^-fluctuations are described by 
the PDF P(x,t\x.i,ti) satisfying equation [j| 

P=-CP with P(x,ti\xi,ti) = 5 XXi , (3.1) 

where C is some linear operator. The principal simpli- 
fication of the problem results from the fact that in the 
Markovian approach (|3.1|) the Q(x,t\xi,U) obeys the SLE 
@: 

Q = - (C + iL)Q , where G(x,t i \x l ,t l ) ^ S XXi , (3.2) 

so that in accordance with eq. (|2.7|) we get for the 
Laplace transform 1Z: 

K= (G) = ((e + iL + C)- 1 ) (3.3) 

The idea of the SLE 1|3.2[1 is based on the simple ob- 
servation that the changes of Q resulted from dynamical 
motion [eq. (|2.2|l ] and stochastic evolution [eq. I|3.1ll ] 
during short time At are written as 

A d G = -(iLg)At and A f G = -{£)g)At, (3.4) 

respectively, so that total change of is given by the 
expression 

Ag = A f g + A d g = -(c + iL))gAt, (3.5) 

equivalent to eq. 13.2fl . 

B. Non-Markovian fluctuations 

1. Continuous time random walk approach 

Non-Markovian y(t)-fluctuations can conveniently be 
described by the CTRWA [which leads to the non- 
Mar kovian SLE 01 f° r G(t)]- It treats fluctuations as a 
sequence of sudden changes of V. The onset of any partic- 
ular change of number j is described by the matrix Pj-\ 
(in {ir}-space) of probabilities not to have any change 
during time t and its derivative Wj—\{t) = —dPj-\{t)/dt 1 
i.e. the PDF-matrix for times of waiting for the change. 
These matrices are diagonal and independent of j at 
J > 1: 

Pj-i(t) = P(t), Wj-i{t) = W(t) = -dP(t)/dt. (3.6) 
For j = 1 

P (t) = Pi(t) and Wo(t) = Wi(t) = -dPi(t)/dt (3.7) 

depend on the problem considered. For non-stationary 
(n) and stationary (s) fluctuations [l5llT6Llr7l| 

Wi(t) = W n (t) = W(t), (3.8) 
Wi(t) = W s (t)=T: 1 J t Q °drW(T), (3.9) 



respectively, where f e = J^dt tW{t) is the matrix of 
average times of waiting for the change 0, 0, . 

It is worth noting some relations for the Laplace trans- 
forms otWj(t) and Pj(t) suitable for our further analysis: 

P,(e) = [1 - WM/e and W s (e) = P(e)/f [Il[ll[l3, 
as well as the representations 

W{e) = [^^(e)]" 1 and P{e) = [e + e/$(e)]- 1 . (3.10) 

in terms of the auxiliary matrix $(e) diagonal in {x}- 
space (see below). 

2. Markovian representation of CTRWA 

In this Section we will briefly discuss the Markovian 
representation of the CTRWA recently proposed in ref. 
[lij . which appears to be very useful for the CTRWA- 
based analysis of the problem under study and, in par- 
ticular, provides the most rigorous method of deriving 
the non-Markovian SLE. 

Suppose that the kinetics of (y — > z/)-transitions in 
{a;}-space is controlled by the Markovian process in an- 
other {<7j}-space, which is governed by the operator A. 
The corresponding PDF <j(j,t) satisfies equation 

a = -Act (3.11) 

describing evolution in {gjj-space and equilibration if the 
operator A has the equilibrium state \e q ) (A|e 9 ) = 0). 
This state is represented as 

K>=E/^'>> (e»l = E 3 -0'l. (WI%) = 1). (3.12) 

where pj are the probabilities of equilibrium population 
of states in {g^j-space. 

The qj -process is assumed to control the evolution in 
{a;}-space as follows: (v — > z/)-transitions occur with the 
rate k„'„ whenever the system visits the transition state 
\t) in {qj}-space. Transitions can lead to the change in 
|j)-state, i.e to (|f) — > \n))— transition with |n) ^ \t). For 
simplicity, we assume that \n) = \t) as well as that A and 
|i)-state are independent of the state in {a;}-space. 

The evolution of the system in {x ® qj}— space is de- 
scribed by the PDF matrix \p) obeying the SLE 

\P) = -{k + iL + k d -K )\p), (3.13) 

where 

K d = k d ® \t){t\ and K = k ® \n){t\ (3.14) 

are the transition matrices (operating in {x <S> qj }-space) 
are diagonal and non-diagonal in {rrj-subspace, respec- 
tively, in which 

K d = y] \v)k, vv (v\, k = y \v)K vv r{v'\. (3.15) 
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and k vv — J2v'hti/) K v'v Equation (|3.13|) should be 
solved with the initial condition 

|p) t=0 = N)E ( 3 - 16 ) 

where \i) = ^jP^H), with (e«K) = Ei4j = L 

The function of interest is the (time) Laplace trans- 
formed PDF in {x}-space 

Q = (e q \p) = X^Pj' in which Pj = O'Ip)- ( 3 - 17 ) 

It is determined by solution |/5(e)} of equation 

\f) = G\i) + GK \$), (3.18) 

where 

G = (fi + A + with VL = e + iL. (3.19) 

The expression for obtained by solving eq. I|3.18|l . 
can be represented in a CTRWA-like form, which in what 
follows [by analogy with the Markovian model (|3.1|l - H3.3|l ] 
will be treated as the solution of the so called non- 
Markovian SLE 0. 

Hereafter, for brevity, we will sometimes omit the ar- 
gument VL of the Laplace transforms of functions under 
study if this does not result in confusions. 

3. Non-Markovian SLE 

Solution of eq. (|3.18(1 leads to the following non- 
Markovian SLE, written in matrix resolvent form |l4j . 

g = Pi(Ci) + p n (A)[i - rWnin^rw^n), (3.20) 

where 

V = k /k d with V V vv > = 1 (3.21) 

* ' V 

is the matrix of jump probabilities with zeroth diagonal 
elements and W ^ (/i = n, i) are written as 

W, L = kddft, with G M = (l + gtkd)- 1 ^, (3.22) 
in which 

^ = (t|(0 + A)- 1 |/i), g t = (t\(Q + A)- 1 \t). (3.23) 

The non-Markovian SLE 13.20|l can also be represented 
in an alternative convenient form |14| 

g = p l + n- 1 $($ + c)- 1 vw l (3.24) 

= P m . +(l- 1 ^ + CY 1 W mz , (3.25) 

where 

t = 1 - V = 1 - k Q /k d (3.26) 



is the Kolmogorov-Feller-type operator describing the 
jump-like motion in {ir}-space, 

* = (.9 f -5n)C 1 + (^3„)- 1 (3.27) 

is the matrix function (diagonal in {x}-space) character- 
izing W(t) [see eq. (l3~TU|) ]: W = W n = (1 + $)-\ and 

W mi = W l /W and P mi = (1 - W mi )/tl (3.28) 

are auxiliary (modified) PDF matrices of the type of wait- 
ing time PDFs. 

The representation (|3.25(l of the non-Markovian SLE is 
somewhat different from that proposed in ref. [3 . Both 
representation are, however, equivalent and can equally 
be used for the analysis, though eq. ( | 3.25|) is closer in its 
form to the CTRWA-like expressions Q. 

The operator C is assumed to have the equilibrium 
eigenstate \e° x ) (i.e. C\e° x ) = 0): 

\4)=J2 P e°»> and {e° x \=Y J <*l- (3-29) 

The eigenstate \e x ), however, is not the true equilibrium 
state of the system under study. The true equilibrium 
state \e x ) is also determined by the behavior of $(e) at 
e — > 0. It is clear from eq. I|3.27|l and definition of the 
PDF matrix Wit) that $(0) = 0, therefore, in general, 
we can write 

$(e) e -° (e/wf, where w = Y, x \x}w x {x\. (3.30) 

In this equation the matrices a (of exponents) and w 
[of characteristic rates (see Sec. V)] are assumed to be 
diagonal in {cc}-space. Moreover, for simplicity, to avoid 
analysis of exotic equilibrium states [l4| we assume that 
a = a is just a parameter independent of x rather than 
matrix. In accordance with formulas l|3.25|l - 13.27|l . for 

such a Q (e + Cw a e 1 ~ a )~ 1 , i.e. the state \e x ) is, 
actually the equilibrium eigenstate of Cw a (Cw a \e x ) — 
0), which is written as 

\e x ) = N- l w~ a \el) with N w - <e°|ur Q |e°). (3.31) 

Possible expressions for \e x ) in some particular models 
of £ are discussed below in Sec. IVA. 

It is important to note that with the eigenstate \e x ) 
the average of any operator Y can be represented as 

(Y) = (e x \Y\e x ). (3.32) 

In particular, as it follows from eq. (|2.ti|l . 

TZ(t) = (e x \g\e x ) = (G). (3.33) 

According to eqs. (|3.20() - (|3.22|l the initial state \i) (in 
{qj}) manifests itself only in the expressions for matrices 
Wi{t) and Pi(t). In particular [HE [HI: 
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a) In the non-stationary n-CTRWA \i) — \n), so that 
Wi = W n , P l = P n and 

Q{tl) = Q n (Q) = A- 1 $(f2)[l>(fi) + C]-\ (3.34) 

b) In the stationary s-CTRWA \i) — |e), where \e) is 
the equilibrium eigenstate [see eq. I|3.12[l ]. consequently 

Wi = W s = P n /r, where f = g n &/pf is the matrix 
of average times (diagonal in {x}-space) with ]f t = (t\e). 
In the considered case \n) = \t), when <J? = l/(kd9t) and 
f = l/(kdPt)- Substitution of W s into eq. H3.20|) yields 



g{Ci) = g s (n) = Ci- 1 - g„(n)£(Of)-\ 



(3.35) 



IV. USEFUL MODELS AND APPROACHES 
A. Models for jump motion 

1. Sudden relaxation model (SRM). 

The SRM Q assumes sudden equilibration in {x}- 
space described by operator 

^=(l-|e )(e |)Q^ 1 , Q = l-T, x P X \x){x\, (4.1) 
in which 

|eo) = E*W, and (e | = (4.2) 

is some auxiliary vector determined by the equilibrium 
vector \e x ): 

\ e x) = J2x p e( x )\ x ) = q\ e o) and (e x \ = (e |, (4.3) 
where 

q = N^Qw- a with N Q = J2 X [(P X - P, 2 )/<]. (4.4) 
According to eqs. I|4.3|l and (|4.4|l : 

\ - y/l-N w-P e (x). (4.5) 



P, 



This relation determines the vector |eo) (in the definition 
of C) which ensures the given equilibrium state |e x ). The 
value of the parameter Nq is fixed by the normalization 
condition for the distribution P x : 



= 1. 



(4.6) 



In the model 1)4. 1|) one gets for any Wi 
K t = (P Qt ) + (q- lp Q )[l - (q-'WQ)}- 1 ^), (4.7) 
where P Qi = (l-W Qi )/Cl, 

W Q = (l + $Q)-\ W Qi =Wi(W Q /W). (4.8) 



The obtained general formulas are simplified in some 
particular models, for example, in the TV-state SRM with 



|e ) = N-^Jx), Q = Q N = l- N-\ 



(4.9) 



In this model q = vT a / Yl T w~ a . It is important to note 
that in the case of two states the model 1)4.9(1 is the only 
possible one. 

Particularly simple expressions are obtained in the 
model of the only characteristic fluctuation rate: 

w = w, for which |e ) = \e x ) = N^Jx). (4.10) 

For s-fluctuations this model assumes the average time 
f e to be independent of x: f e = r e . Together with the 
model (14. 9|) it implies equipopulated equilibration with 
q = 1 so that 

t = Q^Cl-I^Xe^l) and W Q = {1 + Qn^Y 1 . (4.11) 
It predicts, for example, for n-fluctuations (Wi = W) 

K = K n = (P Q )[l-(W Q )}- 1 . (4.12) 

The model l|4.10|l clearly reveals all important specific 
features of the non-Markovian-noise-induced relaxation 
in the simplest form. 

2. Diffusion model 

Another simple model allowing for analytical consid- 
eration of problem is the diffusion model. In this work 
the diffusion model will be applied to the analysis of re- 
combination of the pair of radicals assuming that one 
of radicals undergoes isotropic diffusion in three dimen- 
sional space, say {r}-space, while another radical does 
not move and is located at r = 0. Conventionally, the 
diffusive motion of the moving radical is described by the 
Smoluchowski-type jump operator C = 1 — V |14j |. For 
simplicity we will consider isotropic processes for which 



£ = C D = -A 2 r- 2 V r [r 2 (V r + V r u r )], 



(4.13) 



where r — |r|, V,- = d/dr is the gradient operator, A 2 is 
the average square of the jump length independent of r, 
and u r is the external interaction potential. 

The potential u r is assumed to be of the shape of 
(deep) spherically symmetric square well with the radius 
R much larger than the distance d of closest approach: 
u r = —uq9(R — r) with uq 3> 1. 

It worth noting that within the continuum model of 
stochastic jumps resulting in V(t)- and K(t) -fluctuations 
the corresponding operators V and /C r are just functions 
of r. In the considered spherically symmetric case 



V = 



E \r)V r (r\ and fc = ^ \r)K r {r\ 



(4.14) 
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B. Short correlation time limit (SCTL) 

In practical applications of special importance is the 
SCTL for F(i)-fluctuations in which eq. (|4.7I) can 
markedly be simplified. It corresponds to the large 
characteristic (correlation) rate w c of fluctuations, i.e. 
the large characteristic time of the dependence 3>(Q): 
w c S> ||V||. In our analysis we will discuss the general 
class of models for which &(ft) is represented in the form 
<3E>(f2) = $>((l/w). For such models the correlation rate 
can naturally be defined as w c — \\w\\. 

In the SCTL the relaxation kinetics is described by the 
first terms of the expansion of <f>(fa/w) in small Ct/w c , 

since $(e) is the increasing function of e with $(e) ^5 
[see Sec. IIB and, in particular, eq. (|3.3U|) ]. This fact 
allows one to significantly simplify the problem under 
study. Nevertheless, some important general conclusions 
can be made independently of the form of as it will 
be shown in Sec. V. 



C. Models for quantum evolution. 

Here we describe two systems which will be analyzed 
to illustrate the obtained general results: (1) the isolated 
quantum two- level system (TLS), whose relaxation re- 
sults from V^(i)-fluctuations described within the stochas- 
tic two-state model, and 2) the geminate pair of radicals 
diffusing in the potential well in which spin evolution, af- 
fected by the fluctuating reactivity K(t), in turn strongly 
influences the reactivity. 



c. Simple variants of TLS and two-state model. Be- 
low we will consider two examples of these models: 

1) Diagonal noise 23] : lu s = 0, V x = 0, and V z = 
wq(\x+)(x+\ — therefore 

H^ ± = ±^ (\+)(+\ |— ><—!>; (4.18) 

2) Non-diagonal noise: V z — and V x — v(\x + )(x + \ — 
|ar_)(xZ |), so that 

H V=± =H. ± «(|+)<-| + |-X+|). (4.19) 

In our further analysis the first model is applied to the 
description of dephasing while the second one is used in 
studying population relaxation. 

d. Calculated observables. In the model l|4.16[) de- 
phasing and population relaxation can be characterized 
by two functions: 

1) The spectrum I{ui) which is taken in the form cor- 
responding to Fourier transformed free-induction-decay 
(FTFID) experiments 

7(w) = ^Re{s\n(iu)\s); (4.20) 

2) The difference of level populations 

N{t) = (n\ft(t)\n). (4.21) 
In these two functions 

I s ) = 75ll+->+|-+» and I") = 7= 2 W++)-\—))- ( 4 -22) 

2. Model for reactive radical pairs 



1. Isolated quantum TLS 

a. Two-level model. Quantum evolution of the TLS 
is governed by the hamiltonian (assumed to be real ma- 
trix) with 



H s = \u s a z and V = V z a z + V x a a 



where 



1 

-1 



and er T = 



1 

1 



1+) 



(4.15) 



(4.16) 



b. Two-state model of fluctuations. In the two-state 
model of fluctuations V(t)-modulation is assumed to re- 
sult from jumps between two states (in {x}-space), say, 
\x+) and \x~). It is important to note that in the partic- 
ular case of two states any CTKWA-based kinetics model 
reduces to the simple two-state SRM [see eqs. (|4.9|) and 
H4.10J1 ] in which w = w,Q 2 = 1/2, and 

£ = 2(l-|e a )<e ie |) with | eje ) = i(|a: + ) + |s_». (4.17) 



a. Magnetic field effects. The kinetics of RP recom- 
bination is known to be markedly affected by the RP 
spin evolution which is controlled by the spin Hamilto- 
nian H of the pair. The dependence of the recombina- 
tion kinetics on RP spin state results in a large number 
of phenomena called magnetic field effects [22| . 

In this work we will restrict ourselves to some simple 
and representative effects observed in strong magnetic 
field B for which the Zeeman interaction is much larger 
than the intraradical magnetic interactions (hyperfine in- 
teraction, etc.). We will also consider the effect of the 
external microwave field B\ rotating with the frequency 
lu in the plane perpendicular to the vector B. 

For strong magnetic fields and in the presence of the 
field B\ the spin Hamiltonian governing spin evolution 
of electrons in the pair of radicals, say a and b, can con- 
veniently be written in the frame of reference rotating 
together with B\ with the frequency u> |22^ : 



H 7 — H„ 



H b with = (u li —w)S (1 +tjJiS ll , (4.23) 



where /i = a, b and w M = g^ftB+Y^j Ajlj z is the Zeeman 
frequency of the radical fi possessing some paramagnetic 
nuclei with hyperfine interactions and uj » \{g a + 
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gt,)(3B\. In case of need the matrix representation of 
the Hamiltonian l|4.23[) can be determined either in the 
bases of radical spin states |±) a |±)t or in the basis of 
eigenstates of the total electron spin S = S a + Sj,: singlet 
(\S)) and triplet (\Tq ±)) states, which are expressed as 

\S) = ^(|+>a|->6 - |-)a|+) 6 ), \T ) = ^(|+>a|->6 + 

|-)„|+)6), and \T±) = \±) a \±) b : 

The RP recombination can be treated as a contact re- 
action at a distance of closest approach d using the simple 
model El 



where 



K r = k k s 9(r - d)9(d + A-r), (4.24) 



k s = {P s ,...} (4.25) 



is anticommutator ({V s , p] = V s p + pP s ) in which V s = 
\S) (S\ is the operator of projection on the singlet (\S)) 
spin state of RP. 

The RP is assumed to be initially created in the singlet 
state \S) within the potential well at a distance r = r,i < 
R (fj > d), so that 



p(r,t = 0) = Pi (r) = (4 7 rr l 2 )- 1 (5(r - n)V s 



(4.26) 



b. Observables. In experiments on ma gnet ic field ef- 
fects a number of observables are discussed |22j • Here we 
analyze the most simple ones: 

1) magnetically affected reaction yield (MARY) [2^ . 
measured in the external constant magnetic field, 

2) reaction yield detected magnetic resonance (RY- 
DMR) [22] . i.e. microwave field influenced recombination 
yield. 

In both types of experiments the observables under 
study are recombination (Y r ) and dissociation (Yd) yields: 



Y r = (d/n) 2 AkoTr[P s g(d,r t \e = 0)V S ] (4.27) 

and Yd = 1 — Y r . 

Naturally, the expression l|4.27[) should be averaged 
over nuclear configurations (over lu^). However, in our 
further discussion we will omit this evident procedure 
and analyze the behavior of Y r for fixed oj a and uib (note 
that the case of fixed oj a and tub can be realized exper- 
imentally, for example, with RPs which do not contain 
paramagnetic nuclei). 



V. RESULTS AND DISCUSSION 
A. Isolated quantum TLS 

1. Some general results in the SCTL 

a. Small \\H s \\/w c < 1. Within the SCTL relatively 
simple and general results can be obtained in the case 



\H s \\/w c < 1. In the lowest order in ||$(f2/w c )|| < 1 



K ^TZ n w (u; Q ^ 1 $(fi))/(w Q l>(fi)). (5.1) 

This formula holds for any initial matrix Wi and, in par- 
ticular, for s-fluctuations, if ||f e || ~ l/w c <C 1/||0||. It 
can easily be obtained from eq. (j4.7|l if one takes into 
account that <z -1 Q = N w a . 

b. Large \\H s \\/w c > 1. Somewhat more compli- 
cated SCTL-case ||7J s ||/u; c ~ 1 can be analyzed by ex- 
panding Q in powers of the parameter ^ = ||V||/||iJ s || <^ 
1. In particular, within the general two-level model [eq. 
(|4.1ti|) ] with Vd = in the second order in £ the diagonal 
and non-diagonal elements of p(t) are decoupled and the 
corresponding elements of 7t(t) are expressed in terms of 
the universal function 



V k (t) 

( P \n(t)\p) 



1 

2?ri 



de : : (5.2) 

o e + ke/($(e)) 

"^(0. (A» = «,+- -+),(5-3) 



where 



and 



= ((i\H a \fi), k n = 2Rc(fc+_), 



(5.4) 



= |w7 2 <V„g- 1 [l - W Q (2iu s )]V n ). (5.5) 



2. Anomalous V(t) -fluctuations 

The simplest model for anomalous fluctuations can be 
written as 



$(e) = (e/w) a , (0<a< 1), 



(5.6) 



where w is the matrix of characteristic fluctuation (cor- 
relation) rates diagonal in |x)-basis. For the sake of 
simplicity, w is assumed to be independent of x, i.e. 
W = w(= w c ) so that one can use formula 14.12|) for eval- 
uation of lZ n . The model (|5.6|l describes anomalously 
slow decay of the PDF-matrix W(t) ~ l/t 1+a (very long 
memory effects in the system 9] ) , for which only the case 
of non-stationary (n) fluctuations is physically sensible. 

For small ||iJ s ||/u; c < 1 [see eq. I|5*l|) ] the model l|5.6|) 
yields the expression constituting the important result of 
the work: 



(Cl a - 1 (e))(n a (e)y 1 with n(e) = e + iH. (5.7) 



This expression demonstrates the surprising property of 
relaxation induced by anomalous noise in the SCTL: the 

evolution operator lZ n (e) [and lZ n (t)], which determines 
relaxation kinetics, is independent of the characteristic 
rate w. 
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It is also worth noting that for a — and a = 1 
formula (|5.7I) describes relaxation kinetics corresponding 
to the static and fluctuation narrowing limits pj in which 



where 



(fi-^e)) and TZ n = l/(fi(e)>, 



(5.8) 



respectively. At intermediate values < a < 1 the ki- 
netics is represented by a non-trivial combination of the 
static- and narrowing-like expressions whose relative con- 
tribution is determined by a, as it is seen from eq. <|5.7|1 . 

Of certain interest is the limit a — > 1 in which formula 
(|5.7(l predicts the Bloch-type exponential relaxation: 



n n (e) 



iH B + (a-l){Clhx(Cl))Y 



(5.9) 



The relaxation is controlled by the rate matrix W r = 
(a — l)Rc(f2 ln(f2)), and is accompanied by frequency 
shifts represented by h — i(a~ l)Im(fHn(r2)). The pecu- 
liarity of anomalous relaxation, however, shows itself in 
the independence of matrices W r and h (unlike those in 
the conventional Bloch equation) from the characteristic 
rate w of F(i)-fluctuations. 



3. Anomalous dephasing for diagonal noise 

In the simple model (|4.18(1 of diagonal noise the spec- 
trum can be obtained in analytical form in relatively 
general assumptions on y(i)-fluctuations. 

a. General SRM and SCTL. In the SCTL (i.e. for 
large rate w) within the general SRM (|4.1|) relatively sim- 
ple expression for can be derived without any as- 
sumption on the structure of energy levels (for any num- 
ber of states) 123: 



sin p 



7T (^) 2 + (^) 2 +2^VT COS ^ 



, (5.10) 



where 

p = na and V± = (|w - 2V d fe[±(cj - 2V d )]) (5.11) 

with 6{z) being the Heaviside step-function. In general, 
this formula is too cumbersome for studying the specific 
features of the spectrum. Much more clearly they can be 
revealed with the use of the two-state SRM (see below). 

b. Arbitrary rate w in two-state SRM. The two- 
state SRM <|4.17[) allows for the analytical analysis of the 
spectrum I{u>) for any value of w, i.e. outside the region 
of applicability of the SCTL: 



2* 



ttloq z + z-(n + i]- 1 + 2cos0 o ) 
In this expression 

z± = C(l±w/w ), with e = w(2 1/Q ^), 

n = C( z +)/C( z -) an d O = <p{z + ) + <j>(z-), 



(5.12) 



(5.13) 
(5.14) 



z\ a /y/l + \z\ 2a + 2|z|«cos(<^/2), (5.15) 
sin((p/2) 



CO*) 

4>(z) — sign(z) arctan 



\z\ a + cos(y>/2) 



(5.16) 



Figures 1 and 2 demonstrate a large variety of de- 
pendencies I{ll>) for different values of the parameter 
£ = <jJq/ (2}l a w). In general, the spectrum possesses two 
peaks at to = ±loq which at £ <C 1 and a — * 1 collapse 
into one central peak (see next paragraph). Since in ac- 
cordance with the definition (Sec. IVB) the limit £ -C 1 
corresponds to the SCTL, the collapse of lines for a —> 1 
looks quite natural pj. Simple analysis of eq. (|5.12|) 
and calculation results displayed in Fig. 1 show that in 
this limit general spectrum (|5.12l) reduces to the limit- 
ing one confined in the region |w|/u;o < 1 [see eq. (I5.17f) 
and its discussion]. With the increase of £, however, the 
derealization of the spectrum outside this region is pre- 
dicted (Fig. 2). The derealization becomes especially 
pronounced in the limit a — > 1, as expected. 

c. SCTL within two-state SRM. In the two-state 
SRM (|4.17(l formula (|5.10(l is significantly simplified to 
give the same result as eq. (|5.12|l in the limit 



2 cos if 



(5.17) 



where y = (ojq+oj)/(luq—uj) (see also ref. 11]). Accord- 
ing to this formula anomalous dephasing (unlike conven- 
tional one pj) leads to broadening of I(u>) only in the re- 
gion |w| < ujo and singular behavior of I (to) at ui — > ±luq: 
I {lu) ~ l/(u> ± ajo) 1_Q - For a > a c w 0.59 [a c satisfies 
the relation a c — cos(7ra c /2)] the formula also predicts 
the occurrence of the central peak (at u = 0) of 
Lorenzian shape and width 



wl ~ cos(y>/2)/y a 2 — cos 2 {p/2) 



I{u>) « (27TW0)" 1 tan(v?/2)/[l + (uj/w l ) 2 



(5.18) 



(5.19) 



whose intensity increases with the increase of a — a c (see 
Fig. la). At a ~ 1 the parameters of this peak are 
reproduced by eq. (|5.9|l in which (f21n(fi)) = — (tt/2)lv - 
The origination of the peak indicates the transition from 
the static broadening at a <C 1 to the narrowing one at 
a ~ 1 [see eq. (|5.7f) ]. 

It is worth noting that, of course, for systems with 
complex spectra this transition can strongly be smoothed 
and almost indistinguishable experimentally. 



4- Anomalous dephasing for non-diagonal noise 

a. Dephasing in SCTL for \\H\\/w <C 1. The two- 
state SRM (|4~T7|) enables one to analyze the behav- 
ior of the spectrum L(uj) in the complicated case £ s = 
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lu s /w,v/w <C 1 (within the applicability region for the 
SCTL). In this case dephasing is determined by eq. I|5.7|l . 
After some manipulations one arrives at 



and 



/H = -Re 



1 + r%(iaj) a n_ a (iuj) 



where 



fl (e) = \[{e + 2iE f 



■2iEo)P] 



(5.20) 



(5.21) 



r = 2v/uu s , and E = (w s /2)Jl + r^. (5.22) 



The spectrum I(lu) predicted by eq. (|5.2U|) is depicted 
in Fig. 3 as a function of z = lu/E for two values of 
a and different values of the parameter r. This figure 
demonstrates non-trivial specific features of the shape 
of I (u) depending on the values of a and r. First of 
all, similarly to the case of diagonal noise, (in the limit 
£ s = lu s /w,v/w <C 1) the spectrum I(lu) is localized in 
the region \cu\ < Eq for all values of a. Outside this 
region I(lu) = 0. 

The behavior of I(lu) significantly changes with the 
increase of a: 

1) At small a < 0.6 the spectrum consists of three 
peaks (see Fig. 3a): the central peak at u> = lu c = 
and symmetric edge peaks at lu = lo± = ±Eq (at the 
spectrum edges), in vicinity of which the behavior of 
I(lo) ~ — ^/j| 1_q , where fi — c, ±. The intensity 
of the central peak decreases as a is increased. 

2) For large a > 0.6 each of two symmetric peaks at 
lu = ±Eq split into two ones, so that the spectrum pos- 
sesses five peaks: at ui = u c — 0, at ui = ±E and 
at u) = ±ui a with Lu a < Eq (see Fig. 3b). Moreover, 
the intensity of three original edge and central peaks (at 
lu = zszEq and lu = 0) decreases as a — > 1 while the in- 
tensity two additional peaks increases. Besides, the two 
peaks at lu — ±Lu a approach each other, i.e. Lu a — > with 
the increase of r (as it is shown in Fig. 3b), and in the 
limit r»l collapse into one peak. 

b. Dephasing in SCTL for large lu s /w > 1. The 
model <|4.19[1 reveals some additional specific features of 
the kinetics of phase relaxation in the case of not very 
large w, when £ s = lu s /w > 1. For example, as it is seen 
from eqs. (|5.2|) and (|5.3|l . in the limit ||-ff s || ~ lu s > w 
matrix elements (fi\TZ(t)\fi), (p = H — , — h), which de- 
scribe phase relaxation are given by 



(/W)Im> 



^E a [-k^{wt)% 



where 



E a (-z) 



2ni 



dy— — J 

too y + zy 



(5.23) 



(5.24) 



is the Mittag-Leffler function |26J . In this case the spec- 
trum 

I(lu) — Iq(lu + ) + 7q((jj_), where u>± = lu s ± lu (5.25) 



io(w) 



sin ip 



sin < 



7T \x\(\x\ a + \x\~ a + 2cos0 n ) 

with If = 7TCk, 

X = Lu/(\k + _\ l ' a w), n = {ir\k + _\ l / a w)- 1 , 



and 
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sign(x) arctan 



2 a sin(±^>) 



2 a 1 cos(i^) +LU s /l 



(5.26) 



(5.27) 



(5.28) 



Formulas (|5.25(l - (|5.28|l predict singular behavior of 
I(u>) at lu ~ ±lu s : I(lu) ~ l/|w ± w s | 1_Q , and slow 
decrease of I (to) with the increase of \lu ± lu s \/k 3> 1: 
I(u;) ~ l/|w±w s | 1+a . 

In the limit £ s = lu s /w <C 1 <^(aj) ~ 7ra#(:z:) so that 
7o(w) ~ This means that for £ s <C 1 the spectrum 

I(lu) is localized in the region |a>| < lu s and looks similar 
to I(u) for diagonal dephasing at a < a c and £ <C 1 (see 
Fig. la). For £ s > 1, however, I(lu) is non-zero outside 
this region as well, moreover, in the limit £ s 3> 1 the 
spectrum Iq(lu) becomes symmetric: Iq(lu) — Iq(-lu). 
Such dependence of I(u>) on £ s is also very similar to 
/(w)-dependence on £ found above in the case of diagonal 
noise. 

It is interesting to note that for £ s -C 1 functions 
(/j|7\L(i)|/i) and I(lu) are independent of w [in agreement 
with eq. (|5.7|) ] since fc M ~ (uj s /w) a and k^(wt) a ~ 
(a; s i) a . In the opposite limit, however, fc M ~ w° so that 
the characteristic relaxation time ~ w^ 1 . 



5. Anomalous population relaxation 

Specific features of anomalous population relaxation 
can be analyzed with the model of non-diagonal noise 
(l4~m 

In particular, in the limits \\H S \\ ~ lu s > w and 1 — a <C 
1 with the use of eqs. (|5.2|l . (|5.3|l and (|5.9|l one gets 



N(t) = E a [-k n (wt) a ] and N(t) 



(5.29) 



respectively, where E a (—x) is the Mittag-Leffer function 
defined above and w a = k n (a — > l)w ~ 1 — a. The first of 
these formulas predicts very slow population relaxation 
at t> r r = w- 1 (k n /w) 1/a : N(t) - l/t a . Similar to 7(w) 
the function N(t) is, in fact, independent of w for ^ s = 
lu s /w <C 1 because in this limit k n ~ (dj s /w) Q . In the 
opposite limit > 1 the characteristic time population 
relaxation is ~ w^ 1 since k n is independent of w as in 
the case of phase relaxation. 

For ||7f s ||, ||V|| < w one obtains [H 



i r 



alee e 



+rgr> a _!(e) 
e«+r2fi a (e) ! 



(5.30) 



where £lp(e) and ro are defined in eqs. 15.21|l and 15.22fl . 
respectively. 
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It is easily seen that in the corresponding limits the 
expression (I5.30[l reproduces formulas l|5.29[) with k n ss 
2 a ~ Y cos(ira/2)(E /w) a and w a w tt(1 - a)v 2 /E . Out- 
side these limits N(t) can be evaluated numerically (some 
results are shown in Fig. 4). In general, N(t) is the os- 
cillating function (of frequency ~ Eq) with slowly de- 
creasing average value and oscillation amplitude: for 
t = E t > 1 N(t) ~ l/r Q (except the limit a -> 1). 
At large r = -Eoi one can estimate the asymptotic be- 
havior of N(t): 



where < a < 1, and w and C are the constants with 
( « 1 [small value of C ensures that W(t) > 0]. This 
$(e) corresponds to the waiting time PDF-matrix W(t) 
for which the average time f e = (t) = J^°dttW(t) = w^ 1 
is finite but the higher moments (t") with n > 2 do not 
exist. 

Possible effects of this weakly anomalous noise can be 
analyzed within the SCTL with the use of eqs. (|5.1|) - 
(|5.3() . For example, in the limit ||if || /u> <C 1 one obtains 
formula 



N(r) 



-T(a) 2 Q sin(^)r - 

7T L 



- simp- 



2°^ 



^2 cos(2r + 



(5.31) 



6. Fractional Bloch equation 

The kinetic dependencies found in this section for 
|| V\\/ 1| H s || <C 1 and expressed in terms of the Mittag- 
Leffler function are conveniently represented in the form 
of the equation similar to the conventional Bloch equa- 
tion for the density matrix but with the fractional deriva- 
tives. It is easily seen that the kinetic functions (|5.23() 
and H5.29[) can be considered as solution of the equation 



dH/dt = -iH s n - w a ky^- a Y s ii, 



(5.32) 



where 



r(a) 



lHst Tr / dr . f(r) (5.33) 



is the modified Liouvillc-Ricmann fractional integral op- 
erator, and 



(5.34) 



This equation generalizes the well known classical expres- 
sions [9| to the quantum processes. 



Tl « [e + iH s + Cw- a ((iH) 1+a - (iH s ) 1+a )}- 1 . (5.36) 

which predicts the Bloch-type relaxation of both phase 
and population, but with the rate 



W r = (w- a Re{(iH) L+a - (iH s y +a ) 



(5.37) 



depending on the w as w~ a , i.e. slower than in the con- 
ventional Bloch equation (W r ~ 1/w 0). 

More detailed analysis also demonstrates that in this 
expression for 1Z the terms ~ w(e/w) 1+a are also ex- 
pected. They lead to the inverse power-type asymptotic 
behavior of (/j,\R(t)\fj,) ~ l/t 2+a which, however, is ob- 
served only at very long times t 3> w . This conclusion 
can easily be clarified with the use of the general for- 
mula valid in the case \\H s \\/w < 1, ||^||/||JJ S || < 1 
for any model of $(e). According to this formula for 
the model (|5.35l) the kinetics of both phase and pop- 
ulation relaxation is similar in its mathematical form 
to the probability of fluctuations P(t) [see eq. (|3.6|l ] 
which is, evidently, of power-type behavior at long times: 
P(t) ~ l/t 2+a , in agreement with the above statement. 

Formula (|5.2() allows for making the following general 
conclusion on the relaxation kinetics: the kinetics is of 
anomalous long tailed inverse-power-type behavior for 
any CTRWA-based model of fluctuations assuming sin- 
gular <I>(e) with brunching points. 



B. Anomalous magnetic field effects on RP 
recombination 



7. Weakly anomalous fluctuations 

The analysis presented above shows that the effect of 
the anomalous noise, i.e. fluctuating interaction [whose 
correlation functions is anomalously long tailed: P(t) ~ 
t~ a with < a < 1], on quantum systems can be very 
strong, manifesting itself in anomalous relaxation kinet- 
ics. With the increase of a up to a > 1 the effects of 
anomaly of interaction fluctuations become weaker but, 
nevertheless, they still manifest themselves in the relax- 
ation kinetics. 

To clarify these effects we will briefly discuss the model 
in which 

$(e) = (e/w)+aeM 1+a , (5.35) 



In accordance with eq. I|4.27|l . the observables investi- 
gated in MARY and RYDMR experiments are expressed 
in terms of the operator Q{r, rj|e). Within the CTRWA 
equation for this operator is determined by the model of 
migration of the mobile radical. 

1. Anomalous diffusion and anomalous SLE. 

Here we will consider the anomalous diffusion model 
of migration in which the memory is anomalously long 
and is described by the operator $(e) (I5.t)|) . Notice that 
in the considered process of spatial diffusion the (diago- 
nal) matrix of characteristic rates w in the operator $(£) 
is actually represented as a distance dependent function: 
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w = \r)w r (r\. For simplicity, in our further discussion 
of radical diffusion we assume the rate w r and the expo- 
nent a r to be independent of r: w r = w and a r = a. In 
this case the CTRWA, with $(e)_given by l|5.6[l . is known 
to predict anomalous diffusion |9j. 

The considered anomalous diffusion of the mobile rad- 
ical of the RP, evidently, results in non-stationary fluc- 
tuations of reactivity K r . The effect of these fluctua- 
tion is described by the operator Q(r,ri\e) satisfying the 
non-Markovian SLE (|3.34|) which can be represented in 
terms of the so called fractional diffusion equation (for 
the Laplace transform) as follows 



where 



A 2 (Ad/\ 2 )(k a k s /w) a 
R 2 l + {Ad/X 2 ){k k s /w) a 



> ^ = ^2 e ~"°' ( 5 - 41 ) 



and V s — \S)(S\ is the operator of projection on the 
singlet {\S)) state of the RP. 

Formula H5.40[l is quite suitable for the analysis of the 
problem under study. Unfortunately, in general, in the 
cage model (I5.40|) the expressions for Y r are still fairly 
cumbersome and can be used mainly for numerical esti- 
mations. 



VL r g = -C D {M r g) + 5{r - r % 
where £l r — e + K r + iH and 



M r 



n r $(n r ) = w(fi r /w) 



(5.38) 



(5.39) 



In general, solution of the non-Markovian SLE 
with the Smoluchowski-type Cd [see eq. (I4.13f) ] is a fairly 
complicated (though, in principle, analytically tractable 
[27|]) problem. The most interesting specific features of 
magnetic field effects, however, can quite clearly be il- 
lustrated in some limiting cases allowing for considerable 
simplification of the obtained general expressions. In our 
analysis we will concentrate on one of such cases, corre- 
sponding to the limit of the deep well u r and fast diffusion 
within the this well, which is described by the cage model 

M. 



2. The anomalous-cage model. 

Similar to the case of conventional diffusion, fast 
anomalous diffusion in a large enough and deep potential 
well u r = —uo8(R—r), for which R^> d and uq 3> 1 leads 
to rapid relaxation of the initial non-equilibrium spatial 
distribution of the radical within the well during the time 
td = {R/\) 2 l a /w and formation of the nearly homoge- 
neous quasiequilibrium state (cage) within the well. At 
longer times t > td reaction and dissociation are shown 
to result in the quasistationary decay of this state which 
is, naturally, independent of the distance r% of RP cre- 
ation. At these times, also as for reactions assisted by 
conventional diffusion, the (anomalous) kinetics of the 
process under study is described by lowest pole in the 
expression (I3.34|) (determined by the lowest eigenvalue 
of the operator Cd) H3- In what follows this approxi- 
mation will be called the anomalous-cage model. 

Keeping this lowest pole after some algebraic manipu- 
lations similar to those presented in ref. [27] one arrives 
at formula 



Y r = Tr 



s l r 



1 



l r + l d + {iH z / W y 



3. Specific features of MARY and RYDMR 

To reveal characteristic properties of Independence on 
the parameters of the model we consider simple limiting 
cases which will be somewhat different for MARY and 
RYDMR. 

a. Analysis of MARY In the considered limit of 
strong magnetic field B, in which g^B 3> A^, the ef- 
fect of spin evolution on reaction yield called MARY can 
be studied within the STo-approximation. In this ap- 
proximation, which takes into account that in the strong 
magnetic field limit the contribution of |T±)-terms to the 
reaction yield is negligibly small, the RP spin Hamilto- 
nian is written as 



±6u(\S)(T \ + \T )(S\) with Slu = u a - to b . (5.42) 



(5.40) 



Detailed analysis demonstrates that, in general, the 
anomalous- cage model predicts the MARY-dcpendence 
on parameters of the spin Hamiltonian similar to that 
known in the conventional cage model (conventional 
diffusion assisted processes within the well) [2^ . how- 
ever, with replacement of analytical functions by non- 
analytical ones. 

Most important features of non-analytical MARY- 
dependence on the parameters of the spin Hamiltonian, 
predicted by the anomalous-cage model, can be demon- 
strated in the simple limit of relatively weak magnetic 
interactions: (Suj/w) a <C I4, \\l r \\- In this limit one can 
evaluate MARY with the lowest order of expansion of Y r 
[eq. (|5.40Jl ] in powers of H®: 

Y r w (l s /lo) ~ (l s /l 2 o)Tr[Vs(iH z /w) a V s ], (5.43) 

where l s = Tr(7V r 'P s ) is the reactivity in the singlet 
state and Iq = l s + 14. Straightforward evaluation with 
the use of eq. I|5.43|l gives the expression 

Y r = (ls/lo) - \{l s /ll) cos(7ra/2)|H Q / 2 . (5-44) 

The non-analytical dependence Y r ~ \Suj\ a / 2 is just the 
manifestation of anomalous nature of diffusion within the 
well (anomalous nature of the cage). Noteworthy is that 
in the case a — > 1 the 5w-dependent part vanishes. This 
is because in eq. (|5.43|l the spin dependent contribution 
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to Y r is taken into account in the lowest order in H° 
[~ (H®) a )). At a = 1, however, the term of this order, 
linear in H z , does not contribute to Y r . The non-zero 
contribution, evidently, results only from the second or- 
der term. 

b. Analysis of RYDMR Consideration of the most 
important specific features of RYDMR can significantly 
be simplified in the limit of large STo-coupling u>st = 
5u> ~ (A^) [see eq. Q4.23[l ] and relatively weak microwave 
field ljv {usT /w) a » \\l T \\,l d and < ||Z r ||,/ d . 

In this limit quantum coherence effects on evolution of 
all states with large splitting (~ ^st ) is negligible |27| . 
i.e. their evolution can be treated with balance equations 
(equations for state populations). 

Coherence effects prove to be important only for four 
nearly degenerate pairs of states, of type of the TLS dis- 
cussed above, which describe resonances non-overlapping 
in the considered limit. These four pairs can be combined 
into two groups of pairs of these TLSs: (|±)b|^) a , \T±)) 
and (|±) a |q=)f,, \T±)), denoted hereafter as a± and b±, 
respectively. Transitions in TLS-pairs (fi = a,b), are 
associated with those in corresponding separate radical 
fj,. 

TLS-states |±} a |q=) fc , corresponding to the zeroth z- 
projection of the total spin (S z = 0), are the same for 
systems [i = a and /i = b. However, these systems can 
be considered as uncoupled because in the studied limit 
of large ujst significantly efficient transitions in systems 
a and b occur at different values of uj (i.e. corresponding 
resonances do not overlap, as it was mentioned above). 
For this reason it is possible to distinguish the same states 
(|±) a |=p)f,, belonging to a- and 6-systems, and denote 
them as \a) or |6), respectively (the subscript ± or =p 
can be omitted as it will be explained below). 

The assumed initial population of the singlet state re- 
duces to that of the above-mentioned states \a) and \a) 
(in which S z = 0) with the probability 1/2. Notice that 
these states are reactive. The reactivity matrices are 
similar for all systems and quite accurately determined 
as the two level variant of formula (|5.41l) . These matri- 
ces describe reaction with the same rate approximately 
equal to l s /2, where l s = Tr(V s l r V s ) is the reactivity in 
| iS) -state. 

All the TLSs give the same contribution to the total 
yield Y r , differing only in the resonance frequency (u> a or 
u)b) if they correspond to different radicals. Therefore we 
can combine the identical contributions of the TLSs (J, , 
and \i_ into one Y^ of two times larger magnitude and 
omit subscripts + and — in the notation of the TLSs and 
their parameters, as it has been mentioned above. In so 
doing we arrive at the representation of Y r in the form 



Yr = Y, + Y h 



where Y^, {ji = a, 6), are given by 



(5.45) 



with = and i£ » \l B {V„, ...}. 

For simplicity, we also assume that the microwave field 
uj\ is weak enough so that the effects of the wi-induced 
transitions can be treated perturbatively in the lowest 
order expansion of Y r in \\(H^/w) a /\\l ry d\\ "C 1. 

In these assumptions the yield Y^ can be evaluated 
with approximate expression 

« ±(Z S /Z M ) - \{l s /ll)Tr[V^iHjwYT^ (5.47) 

in which l s = Tr(V s l r 'Ps) and « \l s + Id- 

Calculation using eq. (|5.46|l gives for the magnetic field 
dependent part y r , which is called RYDMR spectrum, 



y r (uj) = y a (u - uj a ) + y b (uj - u b ), 
where y^, (/x = a. b), is written as 

vM = -\{i./iimv ll {%A lx /w) a v^ 

cos(?ra/2) / l s \ 2 (wi/w) Q 
" 2 ^J(l + wVw?)i-«/2- 



(5.48) 

(5.49) 
(5.50) 



It is worth noting some important specific features of 
the RYDMR spectrum (in the case of anomalous diffu- 
sion) predicted by eq. (|5.5U|) : 

1) Unlike conventional Markovian migration, anoma- 
lous diffusion leads to the non-analytical dependence of 
the spectrum on (i.e. on the parameters of spin hamil- 
tonian) which can be obtained in the lowest order in H^. 
Naturally, as in the case of MARY, the lowest-order value 
of RYDMR vanishes in the limit a — > 1 because at a = 1 , 
corresponding to conventional caging, RYDMR ampli- 
tude is determined by the second order term of expansion 
in H^. 

2) At large u> (at line wings) the RYDMR resonance 
contributions y^{uj), (/i = a, b), decrease as y A ,(w) ~ 
1/c<j 2 ~ q , i.e. slower than the Lorenzian line (y(w) ~ 
l/^ 2 ). 

3) The width of resonances in the spectrum is deter- 
mined by the amplitude of microwave field u>\. In other 
words these kind of spectra are always measured in the 
saturation regime . 

4) At first sight, the fact that the width is determined 
by u>i is a consequence of long memory effects on the 
processes governed by anomalous diffusion, i.e the ab- 
sence of the characteristic time in these processes. There- 
fore in the presence of such time caused, for example, by 
the conventional intraradical spin lattice relaxation, the 
width seems to depend on this time. In reality, how- 
ever, this is not true which can easily be demonstrated 
in a simple model assuming this time to result from the 
spin independent decay of radicals with the rate wq. In 
this model the magnetic field dependent yield contribu- 
tions y M , (fx = a, b) are still given by eq. I|5.49|l but with 
(iHn/w) a replaced by [(wq + iH^)/w] a — (wo/w) a : 



= Tr 



W + l d + {iHJw)« 



(5-46) y^) 



1 



Re 



(wo + i^l u 2 + jtf - w% . (5.51) 
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It is easily seen that in the limit wi > too this expres- 
sion reproduces eq. (|5.5U|1 while in the opposite limit 
it predicts ~ 1/(^1 + Lo 2 ) 1 ^ a . Therefore the pres- 

ence of the characteristic relaxation time does not lead to 
the change of the line width of RYDMR spectra slightly 
changing only the line shape. 

VI. CONCLUDING REMARKS 

This work concerns the detailed analysis of the spe- 
cific features of relaxation kinetics in quantum systems 
induced by anomalous noise. Two types of quantum sys- 
tems are considered, as examples: two-level systems and 
radical pairs in a potential well whose recombination is 
assisted by anomalous diffusion. The analysis is made 
with the use of the recently developed convenient and 
powerful method based on the CTRWA and the non- 
Markovian SLE. It demonstrated some important pecu- 
liarities of the kinetics of the processes under study. First 
of all, the relaxation kinetics in both type of systems 



proved to be strongly non-exponential. More subtle pe- 
culiarities were found in spectral characteristics of these 
processes: the line shape, its dependence on the param- 
eters of processes, etc. 

In addition, the non-Markovian SLE proved to be very 
efficient in analyzing not only simple TLSs but also mul- 
tilevel quantum systems. As an example of such systems, 
recombining RP was considered. 

In this work we mainly restricted ourselves to analyt- 
ical analysis of the processes, however, as it was point 
out above, the proposed method also allows for signif- 
icant simplification of numerical treatment of the pro- 
cesses under study, especially in much more complicated 
multilevel quantum systems: magnetic clusters, magnetic 
glasses, etc. 
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Figure captions. 

Fig. 1: The spectrum I{z) = I(lo)uq, where z — lo/luq, 
calculated in the model H4.18fl [using cq. 15.12fl ] for two 
values of £ = uj a /(2 1 / a w): £ = 0.05 (a) and £ = 0.3 (b), 
and different values of a: a — 0.3 (full), (2) a = 0.7 
(dashes), and a = 0.90 (dots). 

Fig. 2: Same as in Fig. 1 but for £ = 0.7 (a) and 
£ = 1.5 (b). 

Fig. 3: The spectrum I(z) = I(u))u>o, where z = uj/Eq, 
calculated in the model l|4.18|l [using eq. I|5.12|l ] for two 
values of a: a — 0.5 (a) and a = 0.9 (b), and different 
values of r = 2v/oj s : r — 0.5 (dots), r — 1.0 (dashes), 
r = 4.0 (full) and r = 8.0 (dash-dots). 

Fig. 4: Population relaxation kinetics N(t), where 
t = Egt, calculated with eq. (|5.30|l (full lines) for two 
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values of ro = 2v/oj s : ro = 1 (a) and r$ = 2 (b), and a = 0.85, and (4) a = 0.75. Straight (dashed) lines 
different values of a: (1) a — 0.95, (2) a — 0.90, (3) represent exponential dependence [eq. (|5.29|) ]. 



